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In our original paper, Ref. we tested the predictions of Shukla et al. in Ref. [2| for 
dense hydrogen, and we came to the following two main conclusions: (1) the prediction 
of a "novel attractive proton-proton potential" is wrong and (2), this error arises from a 
failure of the linearized quantum hydrodynamics (LQHD) model applied by the authors 
outside its applicability range. This conclusion was based on a comparison with our density 
functional (DFT) simulation results. In their comment jsl, Shukla, Eliasson, and Akbari- 
Moghanjoughi (SEA), in fact, have come to very similar conclusions about the validity of 
the LQHD. At the same time there appears to be a misunderstanding about the physical 
basis of DFT that apparantly has led SEA to misinterpret our results. For the benefit of the 
reader we, therefore, give a brief comparison of the two methods below, critically assessing 
their respective strengths and weaknesses. 



To begin with, DFT and the Kohn-Sham equations {4] have been derived rigorously 
from the A^-particle Schrodinger equation of a system of interacting electrons moving in 
the external potential of the nuclei. For practical computations, an approximation has 
to be applied to the exchange-correlation energy functional Excln]- As correctly pointed 
out by SEA, the development of improved functionals Exc[n] is still a matter of active 
current research. The role of E'xcH ^'^ account for the exchange-correlation effects of the 
interacting electrons. All the single-particle quantum mechanical effects, however, are fully 
accounted for by solving the Kohn-Sham equations, which are formally equivalent to a single- 
particle Schrodinger equation for particles moving in an effective potential Vcsi[n];r) that 
is determined self-sonsistently. Thus, as denoted in Table 1 below, DFT correctly captures 
all quantum effects. It becomes accurate at high densities, i.e. in the weak coupling limit, 
when the mean interparticle distance f is much less than the Bohr radius (the scale 
of the local field) sl, i.e. the Brueckner parameter = f/a^ is much less than one. At 
lower densities its accuracy is determined by the exchange-correlation functional Exc- For 
hydrogen reliable expressions for Exc[n] exist which were discussed in Ref. allowing for 
accurate calculations of the proton potential in a dense jellium background. In particular, 
bound states between protons and Friedel oscillations are correctly reproduced. 

Compared to DFT, QHD and LQHD contain three important additional approximations: 
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Property 


DFT 


QHD, LQHD 


(1) 


Quantum diffraction effects 


included for 


approximate, via Bohm potential 






Kohn-Sham states 


("quantum recoil") [3j 


(2) 


Quantum coherence effects 


included for 


missing 




e.g. quantum interference 


Kohn Sham states 


(no phase information) 


(3) 


Spin effects, Pauli principle 


approximate 


approximate, via ideal equation of state 






via Exc 


("Fermi pressure") [3j 


(4) 


Many-particle effects 


approximate 


missing (standard formulation) 






via Exc 


corrections derived from Exc [2, 61 '•'^^ 


(5) 


Accessible temperature 


T = 


T = 


(6) 


Accessible density 


no restriction 


weak coupling, rs ^ 1 








incomplete: misses Friedel oscillations 


(7) 


Resolvable length scales 


no restriction 


I > several f 


(8) 


Strength of perturbation 


no restriction 


weak, (5n < uq 


(9) 


Computational effort 


large 


low, semi-analytical 



TABLE I: Comparison of key properties and limitations of standard density functional theory 
(DFT) and linearized quantum hydrodynamics (LQHD). Exc- exchange-correlation functional, 
Vs = r/aB, is the quantum coupling (Brueckner) parameter, where f denotes the mean interparticle 
distance and ob the Bohr radius, no is the unperturbed density and 5n the induced density 
response. Comments: ^'^^ these references used a simplified version of -Exc from DFT. Extensions 
to finite temperature exist. 

(z) instead of a set of in general complex wave functions, QHD and LQHD solve for real- 
valued hydrodynamic quantities, thereby, losing access to quantum interference effects 
(2), in contrast to the statements of SEA 

(m) as any hydrodynamic theory, QHD averages over a finite volume containing many par- 
ticles, thereby losing the capability to resolve scales on the order of the mean inter- 
particle distance f, cf. point (7), the relevant length scale was determined in Ref. 

{iii) SEA correctly note that QHD is a free electron theory, i.e. it does not contain electron- 
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electron correlation effects and can only be applied to high densities, Vg ^1 (indicated 
by the shaded area in Fig. 1 of Ref. cf. point (6) in the table. Surprisingly, in their 
Letter ^ Shukla and Eliasson extended their claim of the attractive potential to the 
low-density value of ~ 26 supporting it by numerical data. Comparing their data to 
DFT simulations we demonstrated [l| that no "novel" potential minimum other than 
caused by bound states or Friedel oscillations exists. 

In conclusion, even though DFT is - obviously - not an exact theory, as correctly pointed 
out by SEA, it is nonetheless generally more accurate than QHD (with or without lineariza- 
tion) by construction. Therefore, a "failure of DFT" jsl to reproduce predictions of linearized 
quantum hydrodynamics (LQHD) is a serious indication of the failure of the latter. 

This work was supported by the Deutsche Forschungsgemeinschaft via SFB-TR 24 project 
A5. 
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